A&A manuscript no. ?? 

CfPA-95-TH-12 



Your thesaurus codes are: 
12.03.1 12.03.4 



ASTRONOMY 
AND 
ASTROPHYSICS 
14.7.1995 



CMB Anisotropies in the Weak Coupling Limit 

Wayne Hu & Martin White 

Departments of Astronomy and Physics and Center for Particle Astrophysics, 
University of California, Berkeley, CA 94720-7304, USA 

Received ?? ??? 1995/Accepted ?? ??? 1995 



Abstract. We present a new, more powerful and accu- 
rate, analytic treatment of cosmic microwave background 
(CMB) anisotropies in the weakly coupled regime. Three 
applications are presented: gravitational redshifts in a 
time dependent potential, the Doppler effect in reionized 
l/C£narios, and the Vishniac effect. The Vishniac effect can 
Geminate primary anisotropies at small angles even in 
* iat e and minimally reionized models in flat dark-matter 
I flpminated universes with Harrison-Zel'dovich initial con- 
ditions. The techniques developed here refine previous cal- 
culations yielding a larger coherence angle for the Vishniac 
effect and moreover can be applied to non-trivial ioniza- 
,Mi>u histories. These analytic expressions may be used to 
y^jjnpdify results for the standard cold dark matter model to 
fits cosmological constant and reionized extensions without 
Detailed and time consuming recalculation. 
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lG Introduction 

Damping processes in the limit that last scattering of cos- 
mic microwave background (CMB) photons lasted only a 
short duration have been the subject of much recent in- 
terest (see e.g. Hu k Sugiyama 1995a; Bond 1995) and 
require detailed information on the microphysical inter- 
actions that govern photon-baryon coupling (see e.g. Hu, 
Scott, Sugiyama, k Silk 1995). Here we focus on the oppo- 
site limit of a very gradual decoupling of the photons from 
interactions with the matter. This occurs for scattering ef- 
fects in reionized models and gravitational effects between 
last scattering and the present. We significantly improve 
previous analytic treatments of this case (Kaiser 1984, Ef- 
stathiou 1988, Hu k Sugiyama 1994) by calculating the 
anisotropy directly rather than inferring it from the cor- 
responding inhomogeneity of the CMB at last scattering. 

These more accurate calculational techniques are espe- 
cially important in the case of higher order effects which 



can dominate on scales under the thickness of the last 
scattering surface. Here linear contributions have been 
damped away and the anisotropy cannot be obtained by 
the standard procedure of numerically solving the lin- 
earized Boltzmann equation (Wilson k Silk 1982, Bond 
k Efstathiou 1984, Vittorio k Silk 1984). The dominant 
term is the second order Doppler contribution called the 
Vishniac effect (Ostriker k Vishniac 1986; Vishniac 1987). 
Because it is second order, Vishniac contributions are 
strongly weighted toward the present (Hu, Scott, k Silk 
1994) and cannot be considered as a projection of tem- 
perature inhomogeneities at the last scattering surface. 
Indeed, 50% of the signal in temperature fluctuations for a 
reionized cold dark matter (CDM) model comes from red- 
shifts of z < 5. This approaches the redshift at which the 
universe must be reionized to satisfy Gunn-Peterson con- 
straints. Thus even minimally ionized CDM models can be 
expected to have arcminute contributions to temperature 
fluctuations from the Vishniac effect. The standard for- 
malism (Efstathiou k Bond 1987; Efstathiou 1988) cannot 
adequately treat this important case in which last scat- 
tering does not occur on a single well defined surface. We 
evaluate the amplitude and detailed shape of the resul- 
tant minimal Vishniac anisotropy. Although it may well 
be the dominant small scale effect, its presence is unlikely 
to confuse measurements of the sharp damping tail of lin- 
ear contributions. This has important implications for the 
measurement of the geometry and/or thermal history of 
the universe from small scale anisotropies. 



To verify these new analytic techniques, we first de- 
velop them on and apply them to first order contributions 
where comparisons with highly accurate numerical treat- 
ments are possible. We present an analysis of first order 
Doppler effects on small scales which improves upon previ- 
ous analytic work (Kaiser 1984) especially in the late ion- 
ization limit. The weak coupling approximation also accu- 
rately describes even the large angle effects of gravitational 
potential decay in a cosmological constant dominated uni- 
verse. The resultant analytic formulae greatly simplify 
previous treatments (Kofman k Starobinski 1985) without 
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compromising accuracy and allow for a quick and efficient 
calculation of this effect. 



2. First Order Effects 

Let us first verify the weak coupling approximation in 
the case of linear contributions which can be calculated 
precisely by numerical methods. To first order in den- 
sity fluctuations, the Boltzmann equation in Fourier space 
for the evolution of photon temperature perturbations 
AT/T = 0(k, 7, 77) in a flat fio + universe has the 
formal solution 

[0 + *](k,7,r?) =/ o "[(0o + *+7-v 6 )f + 2*] (1) 

[see Hu & Sugiyama 1995b, eqn. (54)] where we have 
neglected the small corrections from the angular depen- 
dence of Thomson scattering and polarization. Here 0o 
is the isotropic temperature fluctuation, \f is the New- 
tonian gravitational potential, v& is the baryon velocity 
with c = 1, k • 7 = kji, and overdots represent deriva- 
tives with respect to conformal time r\ = J dt/a where a 
is the scale factor normalized to the present. The differen- 
tial optical depth to Compton scattering is given explicitly 
by f = x e n e <TTa, where x e is the ionization fraction, n e is 
the electron number density, and a? is the Thomson cross 
section. From left to right, the linear contributions to the 
anisotropy in the integrand are the intrinsic photon tem- 
perature at last scattering (Peebles & Yu 1970), the gravi- 
tational redshift (Sachs & Wolfe 1967), the Doppler effect 
(Sunyaev & Zel'dovich 1970), and the differential gravi- 
tational redshift or integrated Sachs- Wolfe (ISW) effect 
(Sachs & Wolfe 1967). The first three effects contribute 
at last scattering. Consequently, these sources are multi- 
plied by the differential visibility function fe _T which is 
the probability of last scattering in the interval dr\ at r\. 
The ISW effect contributes between last scattering and 
the present as the e~ T suppression at high optical depth 
shows. 

Below the horizon but above the photon diffusion 
scale at last scattering, the photons and baryons evolve 
adiabatically leading to acoustic oscillations as photon 
pressure resists gravitational compression (Peebles & 
Yu 1970). This leads to strong contributions from the 
intrinsic photon temperature and the "Doppler peak" 
structure of anisotropies from standard recombination at 
z* ?s 1000 (see e.g. Bond & Efstathiou 1987; Doroshke- 
vich, Zel'dovich & Sunyaev 1978). However if the uni- 
verse is significantly reionized, last scattering occurs 
much later. In this case, the Compton optical depth r 
drops below unity only after the free electron density 
has been sufficiently decreased by the expansion z* ?a 
100(fi /j7 -25) 1/3 (;c e fi&/i 2 /0.0125) 2 / 3 . At this time, the 
diffusion length approaches the horizon scale. Since diffu- 
sion damps out intrinsic temperature fluctuations in the 



CMB as e _T , acoustic oscillations do not appear in the 
present day anisotropy of strongly reionized models. 

In the reionized case, the other sources in equation (1), 
which arise from fluctuations in the matter, are important. 
These sources can yield strong contributions if last scatter- 
ing occurs after the matter has been released from Comp- 
ton drag. After z<j k, 160(f2o'j 2 ) 1 ^ 5 *e 2 ^ 5 , perturbations 
grow as in linear theory. Here the Hubble constant is given 
by Ho = 100/j km s _1 Mpc -1 . Contributions well after z<j 
can be analytically calculated, following Kaiser (1984), by 
iteratively solving for the feedback of the Doppler term 
into the intrinsic temperature fluctuation at last scatter- 
ing. Averaging equation (1) over angles and keeping the 
leading order term, we obtain for the first order effect 

H 

O + * « / re T iv b j 1 [k(r] - r)')]dr)' « iv b r/k, (2) 



for scales below the horizon at ^, ii; > 1. Here dif- 
fusion and cancellation suppresses the contribution by 
r/k, the optical depth through a wavelength. The £th 
Legendre moment of the temperature fluctuation 0£ = 
i 1 ^ J_ 1 QPi(ji)djj, from the secondary sources is given by 

1 

k 



Qi(k, 770) 



-ii) b T — iv b T + 2k^\ 



(3) 



xe~ T jtikirjo - V )}d V , (£>2) 



which is valid for the Doppler term only under the horizon 
scale at last scattering. 

At sufficiently small scales, the photon traverses many 
wavelengths of the fluctuation during the interval the 
sources in the integral of equation (3) contribute. In this 
case, opposing contributions from crests and troughs tend 
to cancel and damp the anisotropy. Mathematically, this 
can be seen since the spherical Bessel function in equation 
(3) oscillates rapidly in comparison with the source terms. 
Taking the slowly varying quantities out of the integral, 
noting that 

V?r[(£+i)/2] 



ji(x)dx 



(4) 



2 r[(£+2)/2]' 
and expressing the anisotropy in the conventional form of 



9 f°° dh 
Cl =± ^ fc 3 |(MMo)|2 

ft JO K 



we obtain 



C 



(i) 



T[(£+l)/2] 



dk k 3 



(5) 



(6) 



2 \T[(£ + 2)/2]J J k (k Vo ) 6 
x[5 c (77,) + 5i SW (%)] 2 P(fc), 

where the matter power spectrum P(k) = |Aj(fc, ?7o)| 2 and 
the Doppler and ISW sources, 

(7) 



S D ( V ) = (Di + D?U d oe- T , 



5i SW (r?) = 3ff 2 fi - _-- Uge" 



Wayne Hu fe Martin White: CMB Anisotropies in the Weak Coupling Limit 



3 



10- 



10- 



w 10- 

. 1 

^10" 



io- 



io- 




- doppler+prim 

doppler 

- numerical 



0.05 200 
0.1 150 
0.3 100 
0.5 100 



10 



100 



1000 



Fig. 1. Doppler effect in reionized CDM models (Qo = 1, 
h = 0.5, n = 1). For late ionization models 2 r ; < Zd, the an- 
alytic estimate which involves the neglect of Compton drag is 
an excellent approximation to the numerical results at small 
scales. For low flth 2 models, the optical depth between 2 r ; at 
which the universe becomes fully ionized and the present is in- 
sufficient to erase completely the primary signal from standard 
recombination. The total effect is thus described by adding 
e~ 2T of the standard recombination result to the Doppler ef- 
fect. Here and in Figs. 3,4,5 the relative normalization is fixed 
by the amplitude of the matter power spectrum. 



are evaluated at the peak of the £th Bessel function r/i = 
Wo ~ (^+ l/2)/fc. Here we assume linear theory growth for 
the density perturbations, as is appropriate after the drag 
epoch, Aj(fc, 77) = U(77)Aj(fc, 770) an d we have employed 
the continuity and Poisson equations to relate this to the 
velocities and potentials 



fcv 6 (k, 77) 
fc 2 *(fc, 77) 



A 6 k = iDA b (k, 77o)k, 
~(D/a)Hfa A b (k,Tio), 



(8) 



assuming the baryons follow the total matter. Note that 
in an fio = 1 universe D = (77/770) 2 . More generally, it is 
given by D oc H J da/(aH) 3 , where H is the time depen- 
dent Hubble parameter. 

The Doppler and ISW contributions have the same 
angular scaling despite the fact that velocities dominate 
over potentials at small scales (high k) from equation (8) 
(Hu & Sugiyama 1994). This is because the Doppler con- 
tribution suffers severe crest-trough cancellation due to 
the irrotational nature of linear theory flows (Ostriker & 
Vishniac 1986). Crest-trough cancellation does not occur 
if the plane wave oscillates in a direction perpendicular to 
the line of sight k _L 7. However, since vj || k and the 
Doppler effect arises from the velocity component along 



the line of sight, it cannot escape cancellation through 
this mode. Thus ordinarily negligible effects may domi- 
nate small scale anisotropies. 

In Fig. 1, we compare the analytic prediction for the 
Doppler effect with full numerical results. For late ion- 
ized low fi;,/} 2 models, the optical depth between the 
reionization epoch and the present is insufficient to damp 
completely the primary signal from standard recombi- 
nation. To describe the full effect at small angles one 
adds this contribution in quadrature with equation (6): 
£,(tot) _ _|_ e _2r£i(pnm)^ rp n - g s j m pi e expression may 

be useful for post-processing standard recombination re- 
sults to account for reionization without detailed recalcu- 
lation. 
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Fig. 2. Weak coupling approximation for the A ISW effect. 
The analytic approximation based on crest-trough cancellation 
traces the full ISW integral to high accuracy and accounts 
for the falling ISW spectrum. As fio increases, the relative 
contribution of the ISW effect decreases in comparison to the 
Sachs- Wolfe (SW) effect. Here we have chosen a pure power 
law P(k) oc k power spectrum. 



Note that if reionization occurs sufficiently before the 



drag redshift and Slih is low such that 



the 



baryons will not have completely fallen into the cold dark 
matter wells by last scattering. In this case, equation 
(6) overestimates the amplitude but not the angular de- 
pendence of the Doppler effect. Since previous analytic 
treatments {e.g. Kaiser 1984, Hu et al. 1994) also neglect 
Compton drag, we have improved the techniques for all 
cases. 

On the other hand, if fio < 1 the ISW term contributes 
at A-domination, 1 + z\ = (il\/ilo) 1 ^ 3 , and the neglect of 
drag is an excellent approximation for all reasonable fio- 
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Fig. 3. Analytic versus numerical results for A CDM. Compar- 
ison with numerical integration (thick solid lines) shows that at 
large scales, where A ISW and Sachs- Wolfe (SW) contributions 
dominate the anisotropy, analytic formulae (thick dashed lines) 
adequately describe the total anisotropy. At small scales, the 
curves depart from the prediction due to gravitational redshifts 
from the radiation dominated era and acoustic contributions. 
Here we have chosen a high h = 1.0 to minimize these effects. 
Also plotted in thin dotted lines is the analytic prediction for 
a pure P(k) oc k power spectrum with the same normalization 
at large scales. 



Furthermore, since the horizon at A domination is close to 
the present horizon in size, the cancellation approximation 
is valid for all fc's relevant to anisotropies. Thus the above 
approximation provides a simple way of calculating the 
ISW contributions even at large angles for flat A models. 
The combined Sachs- Wolfe and ISW effects can be written 



as Ct = C, 



(SW) 



C. 



(ISW) 



where 



C 



(SW) 



C 



(ISW) 



\ I f k3 {l mk ' v * )l ) 2j2i[HV0 ~ v * )] (9) 

'T[(£ + l)/2]\ 2 f dk k 3 P(k) 



2 \T[(£+2)/2] 



k (kr) ) 



6 '-'ISW 



(m). 



In Fig. 2, we show that equation (9) is an excellent ap- 
proximation of the full integral in equation (3). 

For simple pure power law spectra P(k) = Ak n , the 
Sachs- Wolfe integral can be evaluated analytically and the 
ISW term can be well approximated by a closed form ex- 
pression. In particular, for the n = 1 case the combined 
effect becomes 



£{£ + l)C/ 



2tt 



A 

8^ 



tt4 1.54 



78.5 



r[(i+i)/2] 

m+2)/2] 



(10) 



x(l + 0.08/£) _1 (l 



fi 0.094 )2 .07 



fij- 54 and the 



where we have approximated (fio-D*/a*) 2 
first and second terms in square brackets represent the 
Sachs- Wolfe and ISW effects respectively. For tilts up to 
±0.3 around n = 1, the ISW integral is well fit by mul- 
tiplying the n = 1 result by [(£ + l/2)/£ p ] n ~ 1 , where the 
pivot multipole is £ p = 740% 0A3 h~ 1 . The Sachs- Wolfe re- 
sult can be tilted by employing the standard gamma func- 
tion formula (see e.g. White, Scott, & Silk 1994, eq. 15). 

Realistic models are somewhat more complicated if 
high accuracy is required. In the A CDM model of Fig. 3, 
processing of the initial P(k) oc k spectrum during radi- 
ation domination causes the present power spectrum to 
turn down at small scales. The error incurred by employ- 
ing a power law approximation to the spectrum instead of 
the true processed spectrum is small at the quadrupole, 
but can be significant for somewhat smaller angles. Notice 
that the discrepancy is larger for the ISW effect than the 
Sachs- Wolfe effect since it is generated at later times when 
small scales subtend larger angles on the sky. This clearly 
shows the benefit of calculating anisotropies directly in- 
stead of via the spatial distribution where the differences 
in projection would be lost. We shall encounter this dis- 
tinction again when we calculate the Vishniac effect in the 
next section. 



3. The Vishniac Effect 

Since the first order Doppler effect is mainly cancelled due 
to geometrical reasons associated with irrotational flows, 
second order effects can dominate the anisotropy at small 
angles. While in principle there are many sources in sec- 
ond order perturbation theory, it has been shown that 
the Vishniac effect is the dominant contribution wherever 
such effects are important (Hu et al. 1994, Dodelson & 
Jubas 1995). It is obtained by noting that fluctuations 
in the baryon density Aj change the optical depth r to 
Compton scattering and hence the probability of scat- 
tering. The coupling of the optical depth fluctuation to 
the standard Doppler source yields the Vishniac term: 
7 -q(x) =7 • v 6 (x)A 6 (x). 

Generalizing the Doppler term in equation (1), the for- 
mal solution to the Vishniac effect is 



./o 



(11) 



We have neglected feedback effects from the generation 
of a monopole and quadrupole by the Vishniac effect into 
the full Boltzmann equation. For scales under the width of 
the visibility function fe _T , the Vishniac effect feeds back 
to create a monopole O + \I> = O^qr/k). Just as for the 
first order effect, these contributions are suppressed by the 
optical depth through a wavelength on small scales. For 
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scales much larger than the width, 7 • q sources only a 
dipole. Thus feedback effects are always negligible. 

Following Vishniac (1987), let us decompose the solu- 
tion at the present epoch 770 into spherical harmonics, 



e(k,7,%) = £)e/ ro (k)y /ro (n), 



so that 
|0/™l 2 = 



dilYt m (il) I dr/re T f-qe 




!fc^(jJ-JJo) 



(12) 



.(13) 



Choosing z || k, we note that the azimuthal angle depen- 
dence separates out components of q parallel and perpen- 
dicular to k: 7-q = cos <f> sin 9q± +cos 0q\\ . This dependence 
guarantees that the cross terms between the two compo- 
nents vanish after integrating over azimuthal angles. Thus 
the two contributions add in quadrature and may be con- 
sidered as separate effects. 

Let us consider the parallel component cosflgy = yugy. 
This is to be integrated over the oscillatory function 
e ikii(^-^ ) _ j£ wavelength is much smaller than the 
thickness of the visibility function, the integral will suffer 
severe cancellation unless jj, = 0. Just as in the first order 
case, unless the perturbation is perpendicular to the line of 
sight, we will be looking through many crests and troughs 
of the perturbation. Since the parallel component is pro- 
portional to /i, there is no contribution from this direction, 
leaving a negligible total effect. Because the second order 
parallel term suffers the same cancellation as the first or- 
der term, it can always be neglected in comparison. 

The perpendicular term does not vanish for jj, = and 
thus survives cancellation. Since the final result after sum- 
ming over all k modes has no preferred direction, define 



!©/(*, %)l 2 



4tt 21- 



which corresponds to the former definition of ©£. Using 
the addition theorem for spherical harmonics and the or- 
thogonality of cosm^, we find (Vishniac 1987) 



\&i(k,Vo)\ 2 



1) 



J dnP}[(l-n 



2X11/2 



(15) 



x / drj re T q±e 
'0 



!fc^(jJ-JJo) 



The jj, integral can be performed (Abromowitz & Ste- 
gun 1964, eqs [8.5.1, 10.1.14]), to yield 



|0,(fc,r?o)| 2 = -^+l) 



dllTe q± H^T 



, (16) 



where A 77 = 770 — 77. Notice that this has a simple physical 
interpretation. We know from the spherical decomposition 
that a plane wave perturbation projects onto the shell at 



distance A 77 as ji(kAr]). The projection from k to £ is 
not one-to-one and thus leads to an aliasing which is de- 
scribed by the oscillating tail of the Bessel function. If the 
amplitude of the plane wave is modulated by an angular 
dependence, the projection is modified. In particular, the 
perpendicular component is a face-on projection and suf- 
fers less aliasing. Thus the higher oscillations are damped 
roughly as £/kAr/. 

Now all that remains is to evaluate the q± term. Since 
the Vishniac effect is second order, the convolution theo- 
rem for Fourier transforms tells us that 



q±= (l-£)^E v »( k ')A»(|k-k'| 



(17) 



+v 6 (k - k')A 6 (fc')- 



We can relate the linear theory baryon velocity and den- 
sity by the continuity equation (8). Taking the ensemble 
average of the fluctuation and assuming random phases 
for the underlying linear theory perturbations, we obtain 



<«i(M)«±(M')> = —D{ri)D{ri)D{rf)D{rf) 



(18) 



£p(*')P(|k-k'| 



where P(k) = |A^(fc, /70) | 2 and the projected vector 
k-k' 



I 



kk\ 

~k? ) 



|k-k'| 2 _ 

A bit of straightforward but tedious algebra yields 
1 1 1(1 + 1) 



(19) 



<|0/(*,%)| 2 ) 



(4tt) 2 ril kr] 



-M(k)lf(k)P 2 (k), (20) 



(14) where the mode coupling integral, coming from |d| 2 , is 



M(k) 



dy 



1 dx {l ~ x2){l ~ 2yX)2 

.! (l+t/2-2^) 2 



(21) 



P[k(l + y 2 - 2yxy/ 2 ] P(ky) 



P(k) 



P(k) 



The integration variables arise from y = k' /k and k' • k 
k'kx and the time integral is 



h(k) 



f Vo dn 

/ — Sv^uikAv). 
Jo 



(22) 



Here £V(»7) = DDre~ T t]q/ '(770 — 77), which reduces to 
2(77/77o) 3 fe _T 77o/(l — 77/770) if fio = 1- Note that the inte- 
grand is more strongly weighted toward later times than 
the visibility function itself. It has often been assumed 
(e.g. Efstathiou 1988) that the contributions come mainly 
from last scattering where the optical depth t « 1. This 
causes a significant error in the k to £ projection, i.e. a 
given spatial scale would erroneously be considered to 
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contribute to a smaller angle. It can be a very severe ef- 
fect if the power spectrum peaks toward small scales as in 
the case of baryon isocurvature scenarios (Peebles 1987). 
Even extremely small scales can contribute to observable 
anisotropies at sufficiently late times. 
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Fig. 4. Comparison of different approximations for evaluat- 
ing the Vishniac effect for a standard CDM model (Qo = 1, 
h = 0.5, Q 6 = 0.05) normalized to COBE (Q rms -ps = 20itK). 
The cancellation approximation of equation (23) matches the 
full integral of equation (22) extremely well, whereas the pro- 
jection of the r.m.s. spatial fluctuation (Efstathiou 1988) does 
not. 



The time integral in equation (22) is cumbersome to 
evaluate for high £ due to the spherical Bessel function. 
However just as in the first order case, we can exploit the 
fact that in the cancellation regime the integral extends 
over many wavelengths of the perturbation. Indeed, we 
have just shown that the function Sy multiplying ji is 
even more slowly varying than the wide visibility function 
itself. Taking Sy out of the integral and employing the 
approximation of equation (6), we obtain 



Vi r[(i + i)/2] i 

2)/2] k Vo 
Sv(vi), 




Sy(rit) 



> 1 



(23) 



where recall r/i = r/o — (£ + l/2)/k. This represents an 
excellent approximation for angular scales relevant to the 
Vishniac effect (see Fig. 4). 

The random phase assumption for the underlying lin- 
ear perturbations assures us that there are no cross terms 
between first and second order contributions or different 



k modes. Thus the total anisotropy is obtained by inte- 
grating over all k modes: 



C 



(2) _ 1(1+1) 



(27T)3 „g 



1 r r ]h 

j J T (kn fM(k)ll(k)P\k). (24) 



In Fig. 5, we plot the Vishniac effect for standard CDM 
with late to minimal reionization. Since the Vishniac effect 
is so strongly weighted toward late times, the underlying 
approximation that Compton drag is negligible holds even 
for low fi;,/} 2 models unlike for the Doppler effect. The 
small scale matter power spectrum is well described by an 
analytic fit (Bond k Efstathiou 1984) 



P(k) 



Ak 



[(ak + (bk) 3 / 2 + {ck) 2 Y} 



2/v ' 



(25) 



where a = 6.4(fio/j 2 ) _1 Mpc, b = 3.0(fi /j 2 ) _1 Mpc, 
c = 1.7(f2o^ 2 ) -1 Mpc, and v = 1.13. The normalization 
is related to the CMB ensemble averaged quadrupole by 
A = 7.62 x 10 5 (Q rms _p S /20 y uK) 2 (7i- 1 Mpc) 4 in an fi = 1 
model and by employing equation (10) for a A model. Note 
that C 2 = (47r/5)(Q rms _ PS /2.726K) 2 . 

Most previous calculations (Efstathiou 1988; Hu et 
al. 1994; Chiba, Sugiyama, k Suto 1994; Persi 1995) of the 
Vishniac effect have relied upon the approximations of Ef- 
stathiou k Bond (1987). This involves approximating the 
r.m.s. temperature fluctuation at present, 



l©(Mo)| 



2 

rms 



l P 2 (k) 



*l M (k) r^Ln-Ufsl (26 



and projecting such fluctuations on the sky today via a 
free streaming approximation such as 



£(£■ 



2tt 



2+ 



■k 3 \e(k, Vo )\ 



2 

rms 



(27) 



k=i/( Vo - v ,) 



where 77* is the epoch where optical depth reaches unity. 
This approach is flawed since the fc-space power spectrum 
at the present does not contain sufficient information to re- 
construct the anisotropy in £. As we have seen, anisotropy 
contributions to a given £ are widely distributed along the 
line of sight from last scattering to the present in a form 
that depends sensitively on the power spectrum. In Fig. 4, 
we show that this prior approximation causes a ~ 30% un- 
derestimate of the coherence angle of the Vishniac effect in 
the standard CDM model. The error can be significantly 
more severe in cases where there is more power at small 
scales. Our techniques are also more powerful in that they 
apply to visibility functions with arbitrary features, e.g. 
the bimodal function which occurs in late reionization sce- 



4. Conclusions 

We have presented a formalism for calculating CMB 
anisotropies in the weak coupling or cancellation damped 
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Fig. 5. Primary and Vishniac contriubtions for a range of 
possible ionization histories in the standard CDM model (see 
Fig. 4). Note that even for minimally ionized 2 r ; = 5 — 10, where 
first order anisotropies are nearly indistinguishable from the 
standard recombination case, the Vishniac effect contributes 
a significant fraction of its total in temperature fluctuations. 
Because standard CDM has more small scale power than mea- 
surements suggest, we expect these calculations to be an upper 
limit for CDM-like models. 

regime. When applied to the Doppler effect in reionized 
scenarios, it allows for simple post-processing of the stan- 
dard recombination results. It also gives the A ISW effect 
at large angles to good accuracy. Applied to the Vishniac 
effect, this yields a significantly more accurate description 
of the Vishniac effect than that which has previously ap- 
peared in the literature. As an example of our method we 
have evaluated the contribution from standard CDM for 
a range of reionization histories and shown that the ar- 
cminute scale anisotropies are small but present even in a 
minimally ionized scenario. 

While the small amplitude of arcminute fluctuations 
in standard flat CDM models will make their detection 
difficult in the near future, the fact that they are so small 
implies that their contribution cannot be confused with 
the much larger effect due to the geometry of the uni- 
verse. A positive detection of large CMB fluctuations on 
small scales would be evidence that the geometry of the 
universe is hyperbolic. As we have shown here the sec- 
ond order contributions in models such as CDM (which 
already have more small scale power than observation- 
ally required) are small. Recently, the non-linear analogue 
of the Vishniac effect has been calculated from hydrody- 
namic simulations and shown to be small in comparison 
with the second order effect in CDM (Persi et al. 1995). 
Other higher order contributions have also been calculated 



from simulations. Gravitational lensing only redistributes 
the power, it does not generate more power on small scales 
(Seljak 1995a). It has recently been shown (Seljak 1995b) 
that the Rees-Sciama effect is also small. 

Clusters can also induce anisotropies on the CMB 
from Compton scattering off electrons in the hot cluster 
medium. Hot electrons transfer energy to the microwave 
background leading to spectral distortions in the CMB by 
the Sunyaev-Zel'dovich mechanism. Thus the temperature 
fluctuation will not only have an angular but also a fre- 
quency dependence unlike other sources of anisotropies. 
Ceballos & Barcons (1994) employ an empirically based 
model for clusters. They find that in the Rayleigh-Jeans 
regime, where the Sunyaev-Zel'dovich effect leads to a con- 
stant brightness decrement, the anisotropy at arcminutes 
is on the order AT/T < 10 -7 . Moreover, the signal is in 
large part due to bright and easily identifiable clusters. 
If such known clusters are removed from the sample, the 
anisotropy drops to an entirely negligible level. One there- 
fore expects that the small Vishniac effect dominates the 
small angle anisotropy beyond the damping tail of primary 
anisotropies. Thus the presence of significant arcminute 
scale anisotropies may provide a robust indicator of cur- 
vature in the universe. 
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